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ABSTRACT 

We confirm and characterize the exoplanetary systems Kepler- 445 and Kepler- 446: two mid-M 
dwarf stars, each with multiple, small, short-period transiting planets. Kepler- 445 is a metal-rich 
([Fe/H]=+0.25 ± 0.10) M4 dwarf with three transiting planets, and Kepler- 446 is a metal-poor 
([Fe/H]=-0.30 ± 0.10) M4 dwarf also with three transiting planets. Kepler- 445c is similar to GJ 
1214b: both in planetary radius and the properties of the host star. The Kepler- 446 system is similar 
to the Kepler-42 system: both are metal-poor with large galactic space velocities and three short- 
period, likely-rocky transiting planets that were initially assigned erroneously large planet-to-star 
radius ratios. We independently determined stellar parameters from spectroscopy and searched for 
and fitted the transit light curves for the planets, imposing a strict prior on stellar density in order 
to remove correlations between the fitted impact parameter and planet-to-star radius ratio for short- 
duration transits. Combining Kepler- 445, Kepler- 446 and Kepler-42, and isolating all mid-M dwarf 
stars observed by Kepler with the precision necessary to detect similar systems, we calculate that 
21 U % of mid-M dwarf stars host compact multiples (multiple planets with periods of less than 10 
days) for a wide range of metallicities. We suggest that the inferred planet masses for these systems 
support highly efficient accretion of protoplanetary disk metals by mid-M dwarf protoplanets. 

Subject headings: stars: fundamental parameters — stars: individual (Kepler- 445, Kepler- 446, Kepler- 
42, Barnard’s Star) — stars: late-type — stars: low-mass - stars: planetary systems 


1. INTRODUCTION 

The Kepler-42 exoplanetary system is rather remark¬ 
able: it consists of three sub-Earth-sized planets all or¬ 
biting and transiting a mid-M dwarf host star with pe¬ 
riods of less than two days (jMuirhead et al.l l2012fl , a 
so-called “compact multiple.” Consider that Kepler-42 
is a metal-poor star, wi th a measured [M/H] of -0.27 
in iMuirhead et ail (120121). Fo llowing the calculation of 
iSchlaufman fc Laughlinld2010l l. if we assume the planets 
formed from material in a protoplanetary disk with a to¬ 
tal disk mass equal to 1% of the current host star mass, 
and that the disk had identical metal abundance to the 
star today, that leads to a disk metal content totaling 4.1 
Earth masses. Assuming the three planets have primarily 
rocky compositions, and combining t he predicte d plan¬ 
etary mass-radius relationships of iFortnev et al.l (l2007fl 
with the measured planetary radii in IMuirhead et al.1 
(|2012f) . the total mass of the three planets is calculated 
to be 0.71 Earth masses. Under these admittedly sim¬ 
plistic assumptions, nearly 20% of Kepler- 42’s disk met¬ 
als went into the formation of these three rocky planets. 
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The same calculation for the Sun results in only 5% of 
disk metals contributing to rocky planets (Earth, Venus, 
Mars and Mercury), with significantly more contributing 
to the cores of the solar system’s gas giant planets. 

The preference for metals to contribute to rocky plan¬ 
ets rather than gas-giant cores would be strong evi- 
dence for the planet-formation scenario suggested by 
lLaughlin et all (l2004l l. in which gas-giant-core embryos 
form in the protoplanetary disks around M dwarf stars; 
however, the gas in the disk dissipates before those em¬ 
bryos grow large enough to accrete and are cut-off as 
terrestrial planets. The scenario is already supported 
by the relative scarcity of gas-giant exoplanets found to 
orbit M dwarf stars. Using radial velocity observations, 
iJohnson et al.l (I2010D found a statistical decrease in giant 
planet planet occurrence with decreasing host star mass, 
including M dwarfs in the r adial velocity sample. How¬ 
ever, iCaidos fc Mannl (2014 ) do not find strong support 
for a statistical deficiency of gas giant planets orbiting 
M dwarfs, though they cannot statistically rule out a de¬ 
ficiency. Regardless, the presence of failed embryos in 
some consistent proportion to the amount of available 
metals in the protoplanetary disk would provide support 
for the cut-off accretion scenario. 

It is important to determine whether the planetary 
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system orbiting Kepler-42 is an outlier or a common 
type of system around mid-M dwarf stars, and whether 
metallicity affects the sizes of failed embryos. The occur¬ 
rence of compact multiple systems orbiting mid-M dwarf 
stars has not been thoroughly analyzed. Recent investi¬ 
gations aimed at understanding the pla net occ u rrenc e 
statistics from Kepler, su ch as iHoward et al.l (2012), 
iPetigura et al.l (1201311 and iSilburt et al.l (120 HI) do not 
consid e r mid-M dwa r f star s. iDressing fe Charb onneaul 
(12013D . IGaidos et al.1 (120141 1 and lMorton &; Swiftl (|2014D 
investigated the occurrence of planets around specifically 
M dwarf stars and found a propensity for Earth-sized 
planets, the latter finding roughly two planets per M 
dwarf star with periods of less than 150 days. How¬ 
ever, those investigations did not address the occurrence 
of compact multiples, and many Kepler mid-M dwarf 
stars were not in cl uded in that study, including Kepler- 
42. I Swift et all (I2013H investigated the occurrence of 
compact multiples orbiting Kepler M dwarf stars us¬ 
ing Kepler-32 as a representative of the Kepler M dwarf 
stars. They found that Kepler planet detections around 
M dwarf stars could be recreated reasonably well assum¬ 
ing all planetary systems are clones of the Kepler-32 sys¬ 
tem, indicating a high occurrence of compact multiples. 

In this paper we confirm and characterize two new 
compact multiple systems with mid-M-dwarf host stars: 
Kepler - 4410 and Kepler- 44(0, initially discovered by the 
Kepler exoplanet search pipeline to host two and three 
short-period planets, respectively. Kepler- 445c, Kepler- 
445b, Kepler- 446b and K epler-446 d were first reported 
as planet-candidates by iBurke et al.l (1201411 as KOI- 
2704.01, KOI-2704.02, KOI-2842.01 and KOI-2842.02 re¬ 
spectively, and Ke pler- 446c was repor t ed as a threshold¬ 
crossing event by iTenenbaum et alJ (|2013[ ). and later- 
added to the NASA Exoplanet Archive as KOI-2842.03. 
With regard to the host stars, Kepler -445 was identi¬ 
fied as a mid-M dwarf star i n an o ptical spectro scopic 
survey of late-type KOIs by iMann et al.l (l2013al ). and 
both Kepler -445 and Kepler -446 were identified as mid-M 
dwarf stars in a sepa rate infrared spectros copic survey of 
late-type KOIs bv lMuirhead et aft (|2014fl . We searched 
the light curves for additional planets, and found an addi¬ 
tional planet orbiting Kepler -445 not detected by the Ke¬ 
pler pipeline, which we refer to as Kepler-4456. Refining 
the parameters for the host stars and refitting the planet 
transit light curves, we find that Kepler- 445c is similar to 
GJ 1214b: both are likely mini-Neptunes orbiting metal- 
rich mid-M dwarf stars, and that Kepler -446 is similar 
to Kepler-42 in multiple ways: the low-metallicity of the 
host star, the multiplicity, sizes, and orbital periods of 
the orbiting planets, and the planets’ initial mischarac- 
terization. 

We present the observations, data and analysis of these 
systems in Section [2] In Section [3] we calculate the false¬ 
positive probabilities for the planets orbiting Kepler -445 
and Kepler- 446, confirming the planetary nature of the 
transits. In Section[4]we combine Kepler- 445, Kepler -446 
and Kepler-42 with the full sample of mid-M dwarf stars 
observed by Kepler with similar precision from quarters 
1 to 16 to estimate the occurrence of compact multiple 
systems around mid-M dwarf stars, which we define as 
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two or more planets oribiting with periods of less than 
10 days. Finally, in Section 0 we discuss the implica¬ 
tions of these systems for planet formation scenarios and 
accretion of disk metals by protoplanets orbiting mid-M 
dwarf stars. 


2. DATA AND ANALYSIS 
2.1. High-contrast Imaging 

High-contrast imaging of Kepler planet-candidates 
serves two purposes: with high-contrast images, the level 
of contamination by unresolved objects within the Kepler 
aperture can be determined, and in the event that there 
are no contaminating objects within the aperture, high- 
contrast imaging provides constraints on false-positive 
scen arios that could mimic the detected transit signal 
('e.g.lAd ams et, ahll2012t iLillo Box et al.ll20 12l: lLaw et al.1 

L4U 


l2013t iLillo-Box et al.l 2014 Dressing et al.ll2014f). 


We observed Kepler -445 and Kepler-446 with the 
Keck II telescope at Mauna Kea Observatories using 
the facility near-infrared adaptive-optics imager NIRC2, 
operated in K' band. We observed the stars using 
both conventional imaging and using non-redundant 
aperture masks, which can be placed near an im¬ 
age of the telescope pupil within the NIRC2 cryostat. 
Non-redundant aperture masks enable high-contrast 
imaging by measuring individual spatial frequencies 


of apertures 

Haniff et al.; 1987k (Na 

caiima et al. 1989|: 

Tuthill ct, al. 

1999: Monnier et al. 

20041; Llovd et al. 

2006: Martinache et al. 20091 Kraus & Ireland 20121; 

Ireland 2013|). 

By making the apertures non-redundant, 


high fringe-contrast and phase information is not lost to 
incoherent addition of fringes from redundant baselines, 
enabling contrast performance better than the full-dish 
diffraction limit of the telescope. Phases from combi¬ 
nations of three individual baselines are combined into 
closure-phases, which probe azimuthal asymmetries in 
the image while remaining robust against phase-errors 
from instrument distortions or uncorrected atmospheric 
fluctuations. 

We observed Kepler -445 on UT 17 July 2013 and 
Kepler-446 on UT 29 July 2014. For each target, we ac¬ 
quired three, conventional K' band AO images, each with 
20-second integrations, each dithered across the detector, 
using multiple-correlated sam pling with a Fowler depth 
of 16 dFowler fc Gatlevlll990l f. The non-redundant aper- 
tu re masking image s were acquired in the same manner 
as iKraus fe Ireland! (|2012f ): we used a nine-hole mask, 
with 1.5-meter-equivalent apertures creating baselines 
with separations of 1.5 to 9.2 meters. We acquired six 
aperture-masked images in K band, each with 20-second 
integrations and using mutliple-correlated sampling. 

We reduced the data using the methods described 
in IKraus et al.l (|2008f l. as well as a kernel-pha se tech¬ 
nique with the POISE calibration algorithm (llrelandl 
l2013h . that included subsampling of the Keck pupil 
sub-apertures. The kernel-phase technique yielded 
marginally superior contrast limits, and those limits are 
reported here. 

For Kepler- 445, we detected a faint object at a separa¬ 
tion of 3.908 ± 0.004 arcseconds, with a position angle 
of 279.62 ± 0.05 degrees east of north, and a magnitude 
difference of A I\ l = 7.60 ± 0.07 magnitudes. With such 
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Fig. 1. — Contrast curves for Kepler -445 (top) and Kepler-AAQ 
(bottom) showing the 6cr A K' detection limits for objects at the 
plotted separations. The breaks at 0.150 arcseconds arise from 
a switch from contrast measured using aperture-masking to that 
measured using conventional adaptive-optics imaging and PSF- 
subtraction. We detected a faint object roughly 4 arcseconds from 
Kepler-AA5, with A K = 7.60 ± 0.07 magnitudes. Such a faint 
object could not account for the transits seen around Kepler- 445. 

a large contrast, the object’s contribution to the Kepler 
light curve is negligible and it is unlikely the object show¬ 
ing the transit signals. Nevertheless, we consider the pos¬ 
sibility in our false positive analysis in Section [3j Aside 
from this detected object, we detect no other objects near 
Kepler- 445, with 6er contrast limits shown in Figure [lj 
Interestingly, t his is inconsistent wi t h the Kepler Input 
Catalog (KIC, iBatalha et all l20TTil: iBrown et al.1 1201 If) 
and the NOMAD Catalog (jZacharias et al.ll2005t ). both 
of which record a similar visible-brightness star ~2.”2 
away: KIC 9730159 and NOMAD 1364-0341824, respec¬ 
tively. Considering the infrared adaptive optics imaging, 
and the seeing-limited visible-wavelength imaging we re¬ 
port in the following section which also lacks this nearby 
object, we believe that the NOMAD catalog and the KIC 
may have inadvertently recorded Kepler -445 twice. 

2.2. Host-star Magnitudes 

Both Kepler -445 and Kepler -446 have reliable infrared 
</, H and Ks magnitudes in the Two-Micron All Sky 


Survey ('2MASS. lCutri et~aI1l2003USkrutskie et al.l2006lf . 
They are also included in the KIC, which aimed to mea¬ 
sure SDSS-like u, g , r, i, and z magnitudes for all observ¬ 
able stars in the Kepler field. However, owing to their 
faintness, the KIC lacks u and z magnitudes for Kepler- 
445 and u for Kepler- 446. Other photometric surveys 
of th e Kepler field inclu de the UBV Photometric Sur¬ 
vey (lEverett et all f2012h . which contains both Kepler- 
445 and Kepler -446 and their B and V magnitudes. How¬ 
ever, comparin g the independent V b and measurements 
of Kepler -42 i mMuirhead et ahl (2012:) with the measure¬ 
ments in the UBV Survey shows a significant offset of 
0.23 magnitudes, a 4er discrepancy, which could be due to 
inaccurate calibration effects for very red low-mass stars. 
Yet another survey o f the Keyler f ield is the Kepler- 
INT Survey (KIS. iCreiss et al1l2012D . which contains u, 
g and r magnitudes for Key ler- 446, b ut n o photometry 
for Kepler- 445. Interestingly. ICreiss et al.l (2012 1 found a 
systematic offset of ^0.05 magnitudes between the KIS i 
magnitudes and the magn itudes reported in the KIC, and 
iPinsonneault et ahl (l2012h found a similar offset between 
SDSS and the KIC. We also note that the r’-band Carls- 
berg Meridian Catalogue overlaps significantly with the 
Kepler field, and co ntains a measuremen t for Kepler-AAQ 
( r ’ = 16.731. [Copenhagen University et al.ll2006f) . How¬ 
ever, it does not contain a measurement for Kepler- 445. 

Owing to the discrepancies and incompleteness of the 
various surveys, we chose to measure u , g, r, i and z mag¬ 
nitudes of Kepler -445 and Kepler -446 independently us¬ 
ing the Large Monolithic Imager (LMI) on the 4.3-meter 
Disco very Channel Telesco pe (DCT) in Happy Jack, Ari¬ 
zona dMassev et al.l [20131) . LMI contains the full com¬ 
pliment of SDSS filters and a 6k x 6k full-wafer deep 
depletion e2V CCD, with a field of view of 0.25 x 0.25 
degrees. We observed Kepler-AA5 and Kepler-AA6 on UT 
2014 August 6, along with two nearby SDSS fields cen¬ 
tered on a=290.60299, <5=38.75381 and on a=270.00583, 
<5=270.00583, for calibration. Conditions were photomet¬ 
ric. We observed each SDSS field at four different air- 
masses, spanning 1.00 to 1.30. We observed Kepler-AAb 
and Kepler-AAG twice, at two separate airmasses each, 
with all observations between 1.00 and 1.30 airmasses. 
We also observed Kepler-A2 following the same strategy 
as Kepler-AA6 and Kepler-AA5. 

For each pointing, we cycled through the SDSS filters 
in LMI, acquiring u , g , r, % and z band images using 
2x2 pixel binning. Each exposure was 20 seconds long 
in order to match or exceed the signal-to-noise of the 
archival SDSS photometry without saturating the LMI 
detector. We performed standard flat-fielding and bias 
subtraction, using twilight flats for u and g bands, and 
dome flats for r, i and z bands. To calibrate the LMI 
measured fluxes to SDSS magnitudes, we cross-matched 
all non-saturated point-sources in our SDSS calibratio n 
imag es with the SDSS DR9 catalog of sources (lAhn et al.l 
1201211 . We extracted the flux for each detectable point 
source in the calibration ima ges u sing a modified ver¬ 
sion of DAOPHOT (jStetsonl I1987I 1. and converted the 
object counts into “detector” magnitudes with a arbi¬ 
trary zero point. For those point sources with a reliable 
match in the SDSS DR9 catalog, we calculated the dif¬ 
ference between our detector magnitudes and the SDSS 
archival “PSF” magnitudes, for each band, at each air- 
mass. We only used SDSS PSF magnitudes for objects 
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TABLE 1 

Host-star Magnitudes 


Band 

Kepler-42 

Kepler- 445 

Kepler- 446 

u a 

19.922 ± 0.121 

21.944 ± 0.493 

20.929 ± 0.142 

9 

17.142 ± 0.043 

19.024 ± 0.032 

18.258 ± 0.031 

r 

15.503 ± 0.022 

17.626 ± 0.016 

16.828 ± 0.016 

i 

14.276 ± 0.015 

16.024 ± 0.011 

15.614 ± 0.011 

z 

13.537 ± 0.023 

15.087 ± 0.016 

14.887 ± 0.016 

j b 

12.177 ± 0.021 

13.542 ± 0.029 

13.591 ± 0.021 

H 

11.685 ± 0.018 

12.929 ± 0.035 

13.075 ± 0.026 

K 

11.465 ± 0.018 

12.610 ± 0.028 

12.827 ± 0.024 

Wi u 

11.240 ± 0.023 

12.478 ± 0.024 

12.707 ± 0.023 

W 2 

11.054 ± 0.021 

12.353 ± 0.025 

12.476 ± 0.023 

w 3 

10.831 ± 0.059 

11.252 ± 0.087 

12.931 ± 0.415 


a ugriz measured with DC T-LMI and reported as equivalent 

SDSS PSF AB magnitudes llOke fe Gunn|[l9831 . 

b J, H, and K from 2MASS llCutri et al.lfMCT). 

0 Wi, W 2 and W 3 from WISE jciSn fc et a\\VM3) . 

that were listed as having a star-like point spread func¬ 
tion, not saturated, with a magnitude of less than 22. 
We also limited the calibrating SDSS sources to those 
listed as a star and havin g 0.9 < r — i < 2.25, similar 
to M3 to M6 dwarf stars ilBochanski et al.ll2007f ). Lim¬ 
iting calibration sources to stars of similar spectral type 
eliminates errors arising from different spectral shapes 
across the SDSS bands. This resulted in 6 calibration 
stars for u band, 83 in g band, 249 in r band, 465 in 
i band and 361 in z band for the SDSS calibration field 
with more overall sources. For this field we fitted a line to 
the detector minus SDSS magnitude versus airmass for 
each band. Finally, we applied the fitted relations to the 
measured detector magnitudes for Kepler-AAG, Kepler- 
446 and Kepler M2, again for each band. We estimated 
the uncertainty in our measured magnitudes by combin¬ 
ing the uncertainty due to photon noise with a systematic 
uncertainty based on the linear fit to airmass. The esti¬ 
mated uncertainties are roughly consistent with the dif¬ 
ference between magnitude determinations for the two 
pointings on Kepler-AAG and Kepler-AAG. We averaged 
the two pointings for Kepler-AA5 and Kepler-AAG to ob¬ 
tain our best values, and the results are listed in Table [D 
We do not detect any stellar companions that would sig¬ 
nificantly contaminate the Kepler light curves near either 
Kepler -445 or Kepler -446 in any of the images, despite 
the listed of a similar brightness star near to Kepler -445 
in NOMAD and the KIC. 

With r — i values of 1.68 and 1.28, and i — z val¬ 
ues of 0.86 and 0.64, Kepler -445 and Kepler -446 fall 
squarely in the color-co lor limits for M4 dwarf stars from 
iBochanski et al.l (|2007l Table 1), who compared colors to 
spectral types of M dwarf stars in the SDSS survey. We 
note this only for consistency. In the next section, we 
use our photometry to flux calibrate and stitch together 
optical and infrared spectra, which provide a far greater 
handle on the properties of the host stars. We also note 
that neither Kepler- 42, Kepler -445 nor Kepler -446 shows 
a u-band excess, indicating low chromospheric activity. 
This is further supported by the lack of Ha in emission, 
described in the next section. 

2.3. Moderate-resolution Spectra and Stellar Parameters 

Neither Kepler -445 nor Kepler -446 currently has a pub¬ 
lished astrometric parallax measurement. With astro¬ 


metric parallaxes, relatively accurate mass-luminosity re¬ 
lations can be used to determine ste llar masses b ased on 
their absolute infrared ma gnitudes (|Henrv fe McCarthy! 
119931 : iDelfosse et al.l 120061 1 . The stellar masses can be 
combined with empirical or theoretical mas s-radius rela¬ 
tions hips to determine stellar radius (e.g. iTorres et all 
l20lQlf . Or, in some cases, stellar parameters can be 
determined from the transit light curve itself, either 
by measuring photo metric signatures from asteroseis¬ 
mic pulsations (e.g. iBedding et al.l 120111 : iHuber et al.l 
SoH), or by using an exoplanet transit light curve to 
infer the host star’s den sity assuming a low-eccentricity 
or circular orbit (e.g. iSeager fe Mallen-Ornelasl l2003t 
iCarter et al.ll201lH or a combination of both techniques 
Ve.g. iBallard et aill2014H . 

In the case of Kepler -445 and Kepler- 446, asteroseis- 
mic signatures are very difficult to measure because, be¬ 
ing M dwarf stars, they are dense, with mean densities 
over 10 gm/c, where the asteroseismic signals are low 
and the stellar oscillation frequencies are high. Also, 
the signal-to-noise of the Kepler light curves is simply 
not high enough to provide strong constraints on the 
stellar density from transit light curve fitting. With¬ 
out astrometric parallaxes, asteroseismic or transit light 
curve constraints, we must rely on colors and spec¬ 
troscopy to determine the stellar parameters. In the 
case of M dwarf stars burning hydrogen on the main- 
sequence, spectroscopy probes primarily effective tem¬ 
perature, T e g , and metallicity, [Fe/H] or [M/H], Typi¬ 
cally, surface gravity, log(g), is another physical parame¬ 
ter probed by spectroscopy. In the case of M dwarf stars 
on the main sequence, log(g) is predicted to be a strict 
function of T e ff and metallicity. With T e ff and metallicity 
alone, the stellar mass, radius and bolometric luminosity 
can be determined using either predic tions from stellar 
evolutionary mod els fe.g. iDotter et all l2008li . or empir¬ 
ical r elations (e.g. iSegransan et al.ll2003t iBovaiian et alJ 

IMD. 

Moderate-resolution infrared spectra (R~2700, 1.5 to 
2.5 nm) for Kep l er -445 and Kepler -446 were obtained by 
iMuirhead et all (|2014T ) . who determi ned stellar parame¬ 
ters b ased on the K-band indices of iRoias-Avala et al.l 
(|2012f) . then interpolated those values on a new set 
of models based on the Dart mouth Stellar Evolution 
Database (|Dotter et all I2008H . calculated by G. Fei- 
den. A visible low-resolution spectrum for Kepler -445 
R~900, 3200 to 9200 A) was obtain by iMann et al.l 
2013aD , who used a new calibration to measure ef¬ 
fective temperature. T hey then use d the mass-radius- 
temperature relations of IBovaiian et all ( 2012 ) for their 
sample, but were unable to determine the stellar mass 
and radius of Kepler -445 as it was too cool for these re¬ 
lations. Kepler -446 was not discovered by the Kepler 
pipeline at the time of their study. 

In this paper, we c ombine the te c hniqu es from 
IMann et all (l2013al ) and IMuirhead et all (|2014l ) to pro¬ 
duce the best possible phy sical para meters of the stars. 
We choose to use the IMann et all (l2013a i visible-light 
spectroscopic effective temperature calibration, since it is 
calibrated using truly empirical measurements of nearby 
M dwarf stars from o ptical-long baseline interferometry 
(|Bovaiian et al.1120121 ). We use t he K-band spectroscopi c 
metallicity measurements from IMuirhead et alJ (|2014fl . 
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as tho se used empirical calibrations bv lRoias-Avala et al.1 
(12012D . who verified the metallicities by comparison to 
space-motions of nearby M dwarf stars. Finally, we inter¬ 
polate those values onto the new suite of Dartmouth evo - 
lutionary models as was done in iMuirhead et all . ( 201 41. 
but with a small correction based on the empirical mea¬ 
surem ents of Barnard’s Star (similar to IMuirhea d et al.l 

[20l2ll . 

To do this, we used primarily the archival spectroscopy 
described above, but acquired a new visible spectrum of 
Kepler- 446. Our additional spectrum of Kepler -446 was 
acquired with the SNIFS instrument on the University 
of Hawaii 88-inch Telescope at Mauna Kea Observatory. 
We followed identica l observing and reduction procedures 
described in lMann et al.1 (|2013aT) . 

To stitch the visible and near-infrared spectra together 
for each star, we combine the r magnitudes measured 
in the previous section with the H and Ks magnitudes 
from 2MASS. We acquired the filter transmission curves 
for each band (r, H and K) and computed a synethic 
spectral magnitude using the transmission curve and the 
measured spectrum, keeping in mind that 2MASS used 
a modified Vega system, and our measured_SDSS magni¬ 
tudes use the AB system (Oke fc On mil! lbs: };.. We calcu¬ 
lated the difference between the synthetic spectral mag¬ 
nitude and the measured magnitude from each respective 
survey, and used that difference to renormalize and com¬ 
bine the visible and near-infrared spectra. 

Figure [2] plots the visible and infrared spectra of 
Aep/er-445 and Kepler- 446, with spectra of Barnard’s 
Star from IMann et al.1 (|2013al) and Kepler -42 from 
IMuirhead et al.l ( 20141) for comparison, subject to an ar¬ 
bitrarily total normalization such that they overlap in H- 
band. Clearly, the spectra are very similar: are all identi- 
fiably M4 dwarf stars. We did not correct for interstellar 
reddening because all stars are likely well within 150 pc 
of the sun (see Section [2]) . None of the stars show Ha in 
emission, indicating low quiescent-activit y and tha;Uall 
of the stars are likely old (t > 5 Gvr. lWest et al.ll2008D . 
However, we note that Barnard’s S tar d oes show occa¬ 
sional Ha flaring fe.g. lPaulson et al]l2006l ). and that age- 
activity relationships are a statistical tool for estimating 
the age of an ensemble of stars, and are subject to large 
uncertainties when applied to any specific star. 

To determine any offsets to apply to the new Dart¬ 
mouth models, we interpolated the empirical effective 
temperature and metallicity measurements of Barnard’s 
Star onto the new Dartmouth models to determine the 
predicted stellar mass and radius. The predicted mass 
and radius are marginally different from the empirical 
mass (determined u sing mass-luminosity relationships of 
Dclfoss e et al.ll200nh and the empirical radius (measured 
using optical long-baselin e interferometry, lLane et ahl 
l2001t iBovaiian et al.ll2012l ). with offsets of -0.004 Msun 
and 0.002 Rsun , both of which are well below the typi¬ 
cal uncertainty in those quantities. Nevertheless we still 
applied those corrections to the interpolated values. 

Our results for the stellar parameters appear in Table 
[2] In addition to calculating the parameters for Kepler- 
445 and Kepler- 446, we also recalculated the parameters 
for Kepler- 42 us ing i dentic al met hods, and find agree¬ 
ment with IMuirhead et al.1 (I2012T ). who used a similar 
technique, but applie d a correcti on to the original Dart¬ 
mouth models dDotter et al.l l2008h . rather than the new 


set used in IMuirhead et al.1 (120141) and this study. We es¬ 
timated the distance to the stars by invert ing the mass- 
luminosity relations of lDelfosse et al.l (120001) to determine 
the stars absolute A'-band magnitudes, and compared 
that to the measured A'-band magnitude from 2MASS. 
All three stars are relatively nearby with distances less 
than 150 parsecs. 

The stellar parameter determinations for Kepler-445 
and Kepler-446 presented here revise values publishe d 
in the literature. Compared to IMuirhead et all ( 2014 1, 
the masses and radii presented here a re slightly larger 
owing to the use of the IMann et~aD (j2013a ) effective 
temperature calibration rather than the K-band tech¬ 
nique. The most recent evaluation of stellar param¬ 
eters for Kepler targets was compiled by iHuber et al.1 
([20Ti . who assign both Kepler -445 and Kepler -446 
smaller and larger radii (respectively) than presented 
here. We ascribe this to erroneous metallicity deter¬ 
minations in that catalog, which lists Kepler-445 as a 
metal-poor star ([Fe/H] = -0.380) and Kepler-446 as 
a metal-rich star ([Fe/H]=+0.30), the inverse of our 
measur ements. The sourc e for the Kepler-446 parame¬ 
ters in IHuber et al.l ( 2014 ) is iDressing fc Charbonneaui 
(|2013l) . w ho used photometry to estimate s tellar param¬ 
eters. As IDressing fc Charbonneaui (120131) admit, pho¬ 
tometry alone places poor constraints on M dwarf metal- 
licity, and Kepler-446's erroneou sly l arge radius in both 
IDressing fc Charbonneaui (I2013D and IHuber et al.1 (12013D 
is a direct consequence of the imprecise metallicity mea¬ 
surement. 

2.4. High-resolution Spectra and Space Motions 

We also acquired higher-resolution spectra in order to 
measure the absolute radial velocities of the stars and 
their galactic space motions. We observed Kepler-445 
and Kepler-446 with the Echellette S pectrogra ph and Im¬ 
ager (ESI) on the Keck-II telescope (Shci nis et, al.l120021) 
on UT 25 July 2014. We operated ESI in echellette 
mode using a slit width of 0.5 arcseconds, achieving a 
resolving power of ~ 8000 from 4000 to 10000 A, spread 
across ten cross-dispersed orders. For Kepler-446 we ac¬ 
quired a single 1000 second exposure, and for Kepler -445 
we co-added three 900 second exposures. We also ob¬ 
served GJ 687 and GJ 905 as M-dwarf radial velocity 
standard sta rs with published ra dial velocities and un¬ 
certai nties in lNidever et al.l ()2002l ) and iDeshoande et al.1 
(2012), respectively. We reduced the data using the 
publicly available ESIRedux pipeline^, using dome flats 
for flat-fielding and arc lamps for wavelength calibration 
(|Prochaska et alJf2003t iBochanski et alJl2009ft . 

We then cross-correlated the radial velocity standard 
spectra with Kepler-445 and Kepler-446 to measure their 
absolute radial velocities. We used the full-width-at-half- 
maximum of the peak in the cross-correlation function as 
our measurement uncertainty, combined in quadrature 
with the archival measurement errors for the radial ve¬ 
locity standard stars. We found agreement within our 
uncertainties when using GJ 687 or GJ 905 separately 
as standards, and we report the mean of those measure¬ 
ments here. For Kepler-445 we measured an absolute 
radial velocity of -61 ± 1 km s -1 , and for Kepler-446 we 
measured -118 ± 1 km s -1 . 

5 http://www2.keck.hawaii.edu/inst/esi/ESIRedux/index.html 
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Fig. 2. — Visible and near-infrared spectra of Kepler-446 and Kepler- 446, with spectra of Barnard’s Star and Kepler-42 for comparison, 
all normalized to match in H-ba nd (A ~ 1.65 pm). The visible spectrum of K epler- 446 is part of this work, with the remaining visible 
and infrared spectra coming from 1 Mann et a l. (2013a) and Muirhead ct al. (2014 ). Th e spectra were first ren ormalized piece-wise to match 
infrared photometry from the Two-Micron All Sky Survey f2MASS, ICutri et al . 2003; S krutskie et al.l 120061) and our independent SDSS r 
band photometry. The similar shapes to the spectra strongly support our interpretation that Kepler- 445 and Kepler-446 are mid-M dwarf 
stars with similar masses and radii to Barnard’s Star and Kepler-42. We have made no correction for interstellar reddening, as these stars 
are all expected to be within 150 pc of the sun (see Table 


TABLE 2 

Stellar Parameters 1 " 1 


Mid-M Dwarf Star 

TeS 

[Fe/H] 

[M/H] 

M* 

R* 

Distance 

Kepler- 445 
Kepler-446 
Kepler-42 

Barnard’s Star 

3157 ± 60 K a 
3359 ± 60 K 
3241 ± 57 K a 
3238 ± 11 K a 

+0.27 ± 0.13 c 
-0.30 ± 0.12 c 
-0.48 ± 0.12 d 
-0.39 ± 0.17 c 

+0.19 ± 0.12 c 
-0.21 ± 0.12 c 
-0.33 ± 0.12 d 
-0.27 ± 0.12 c 

0.18 ± 0.04 Mq 
0.22 ± 0.05 Mq 
0.15 ± 0.03 M 0 d 
0.159 ± 0.013 M 0 a 

0.21 ± 0.03 Rq 

0.24 ± 0.04 Rq 

0.18 ± 0.02 R e d 
0.1867 ± 0.0012 R 0 e 

~ 90 pc 
~ 120 pc 
~ 40 pc d 

1.824 ± 0.005 pc f 


Va lues w ithout a mark ar e from this work, values with a mar k are from the followi n g: (a) IMann et al.l J20 1 3al ) , Ibl [Rojas- Ayal a, et alj 
( |20 1 211 . {c) IMuirhead et al.1 120 1 41 ) . (d) IMuirhead et alj (120 1 21 1. (e) IBovaiian et alJ 120 121 0 (f) Hipparcos astrometric parallax from 
Ivan Leeuwenl <20071 1 ■ 


When combined with archival proper motion measure¬ 
ments and our distance estimates from the previous sec¬ 
tion, we can estimate the stars’ space motions. The 
most recent proper motion measurements available for 
Kepler^ 445_aud Kepler - 446 are from the PPMXL cata¬ 
log dRoeser et al.l l201Ch . They report a proper motion 
of /i a =42.2 ns= 132.7 mas yr -1 for Kepler-M5 and /i a =- 
13.2 ns=- 30.6 mas yr -1 for Kepler- 446. Both proper mo¬ 
tion measurements are relatively large and are consistent 
with M dwarf stars within 150 pc. 

Combining proper motion, radial velocity and distance 
estimates, we measure the galactic space velocity mo¬ 
tions of Kepler -445 to be U = 59, V = -39, W = 9 km 
s -1 , and for Kepler -446 we measure U = 8, V = -98, 
W = -31 km s” 1 . We c orrected for the galactic solar 
motion using the values oflCoskuno glu et all (l2011h :U = 
-8.5, V = 13.38, and W = 6.49 km/s. Figure [3] plots the 
space motions for Kepler- 445, Kepler- 446, Kepler-42 and 
Barnard’s Star with nearby dwarf stars in a Toomre dia¬ 
gram w ith thin disk, thick disk and halo boundaries from 
iFulirmaiml (120041) . Kepler -42 and Kepler -446 are located 
in the thick disk regime, consistent with their low metal- 
licities, whereas Kepler-445 is located in the thin disk 
regime, consistent with its high metallicity. We note that 
the distances to Kepler-42, Kepler-445 and Kepler -446 
and their corresponding tangential motions are highly 
uncertain, as is the true solar motion through the galaxy. 
Therefore, we estimate the uncertainty in their galactic 


space motions to be ~10 km s . 

2.5. Kepler Photometry 

Kepler-445 was observed by Kepler in quarters 6, 8 and 
9 in long-cadence mode as part of guest observer pro¬ 
gram G020031, a search for microlensing by observing 
high-proper motion stars dDi Stefanoll201Cl) . It was later 
observed in quarters 12, 13, 14 and 16 in long-cadence 
mode and in quarter 17 in short-cadence mode as an 
exoplanet search target as part of the primary Kepler 
Mission. However, quarter 17 was cut short due to a 
malfunction in the spacecraft ending the primary mis¬ 
sion, so we did not include it in our analysis. 

Kepler -446 was observed by Kepler in quarter 7 in 
long-cadence mode as part of guest observer program 
G020001, a search for eclipsing binar y systems in a sam¬ 
ple of M dwarf stars dHarrisonl 120101) . It was again ob¬ 
served in quarters 12, 13, 14 15 and 16 in long-cadence 
mode and in quarter 17 in short-cadence mode as an 
exoplanet search target as part of the primary Kepler 
Mission, though we do not include the quarter 17 data 
for the reason stated above. 

2.5.1. Known Transiting Planets 

At the time this manuscri pt was su bm itted, the NASA 
Exoplanet Archive (NEA, lAkeson et alJl20l3l) listed two 
planets orbiting Kepler-445 (KOI-2704), with orbital pe¬ 
riods (P) of 4.871229 ± l.le-05 and 2.984151 ± l.le- 
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Fig. 3.— Toomre Diagram of G K and M type stars with mea¬ 
sured trigonometric parallaxes greater than 100 mas , the same 
stars used in a similar figure in [Muirhead et aLl (120121 , , Fig. 9). 
We includ e the thin disk, thick disk and halo boundaries from 
IFuhrmannl (120040 for reference. Barnard’s Star, Kepler- 42, Kepler- 
445 and Kepler- 446 are all shown. Kepler-42 and Kepler- 446 are 
located in the thick disk regime, consistent with their low metal- 
licities, whereas Kepler- 445 is located in the thin disk regime, con¬ 
sistent with its high metallicit y. We correcte d for the solar motion 
using the values of ICoskunoglu et al.l (|201lQ :U = -8.5, V = 13.38, 
and W = 6.49 km/s. We note that the distances to Kepler-42, 
Kepler-445 and Kepler- 446 and their corresponding tangential mo¬ 
tions are highly uncertain, as is the true solar motion through the 
galaxy. 
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05 days, for Kepler-AA5c and Kepler- 445b respectively 
(KOI-2704.01 and KOI-2704.02), with planetary radii 
(Rp) of 2.9 and 1.96 i?®, respectively. The NEA listed 
three planets orbiting Kepler -446 (KOI-2842) with or¬ 
bital periods of 1.5654088 ± 3.3e-06, 5.148921 ± 2.2e-05, 
3.03617925 ± 5.49e-06 days for ifepfer-446b, Kepler- 446d 
and Kepler- 446c respectively (KOI-2842.01, KOI-2842.02 
and KOI-2842.03), with planetary radii (Rp) of 25 ± 
15, 26 ± 15 and 2.393 ± 0.582 R e . For Kepler- 446, the 
NEA listed notably large impact parameters ( b ) of of 1.17 
± 0.92, 1.19 ± 0.99 and 0.90 ± 0.177, for Kepler- 446b, 
Kepler- 446d and Kepler- 446c respectively. 

The planetary radii listed in the NEA as initially 
downloaded were determined using either an MCMC 
fitting routine, as in the case of A epfer-445c, Keyler - 
445b, Kepler-AA&o and Kepler- 446d (jBurke et al.l 12014 ). 
or by the Kepler pipeline when initially searching for 
threshold-crossing events , as in the case of Kepler- 446c 
dTenenba um et ahl 120141. with the stellar parameters 
from iHuber et all (|2014f) . For Kepler- 446b and Kepler- 
446d We attribute the large radii and uncertainties in the 
NEA to correlations between the transit impact param¬ 
eter and planet-to-star radius ratio. The parameters are 
typically not strongly correlated; however in the case of 
Kepler -446 the transit durations (T) are relatively short, 
lasting only about an hour each, and the long-cadence 
integrations times are ~30 minutes. This results in the 
light-curve being significantly smoothed and appearing 
to be V-shaped, where a truly flat-bottomed, low-impact 
parameter transit curve is indistinguishable from a graz¬ 
ing eclipse. 

We note that when this manuscript was accepted, the 
parameters listed in the NEA had changed, listing more 
reasonable values for the planet radii for Kepler- 446, but 


still significantly larger than the values we calculate in 
Section 12.5.31 

2.5.2. Independent Planet Search 

The Kepler -445 and Kepler -446 transiting planet- 
candidates found by the Kepler pipeline did not uti¬ 
lize all quarters of data available currently. Therefore, 
we conducted our own search for additional transiting 
planets on the complete, available Kepler dataset for 
Kepler- 445, Kepler-AA6 and Kepler-A2. We downloaded 
the Kepler light curves for Kepler -445 and Kepler-AAQ 
from the NASA’s Mikulski Archive for Space Telescopes 
(MAST). The light curve data files contain flux values 
measured using simple aperture photometry on the Ke¬ 
pler pixels using pre-defined apertures (SAP_FLUX), and 
flux values that have been detrended using custom tools 
to remove slowly-varying instrumental fluctuations in the 
Kepler photometric response as well as fluctuations due 
to stellar rotation and intrinsic variability: the “pre¬ 
search data conditi oning simple aperture photometry 
flux” ('PDCSAP_FLUX. iSmith et ahll2012l) . We choose to 
use the PDCSAP_FLUX for the independent planet search 
and for fitting the transit events. 

For each star, we fitted the PDCSAP_FLUX light curve 
with a cubic basis spline (cubic B-spline) with break¬ 
points located 1.5 days apart, and divided the light curve 
by the B-spline fit to remove stellar and instrumental 
variability. We excluded outlier data points (both astro- 
physical and otherwise) from the B-spline fit by itera¬ 
tively fitting the B-spline, locating points falling more 
than 3cr away from the fit, and re-calculating the B- 
spline while excluding the outliers. We repeated this 
process until convergence, typically 5 iterations. We 
then searched for transits by calculatin g a Box Least 
Squared periodogram (BLS, iKovacs et alJl2002l ) for each 
star. We evaluated the BLS periodogram over periods 
ranging from 0.15 days (or 3.6 hours) to the total du¬ 
ration of the Kepler observations. We evaluated the 
BLS power spectrum at roughly 10 5 to 10 6 discrete pe¬ 
riods, depending on the total time baseline of Kepler 
observations, and we spaced the trial periods to ensure 
that the BLS signal of a short duration transit around 
a mid M-dwarf would not be smeared out by coarse pe¬ 
riod spacing. After we calculated the BLS power spec¬ 
trum, we subtracted away a noise floor from the power 
spectrum and estimated its typical scatter by calculating 
the Median Absolute Deviation (MAD) and dividing by 
0.67 to convert to an equivalent standard deviation. We 
considered any peak with a BLS signal-to-noise ratio of 
greater than 9 to be significant. Upon detecting a sig¬ 
nificant signal, we masked out the signal in question and 
re-calculated the periodogram. 

In all three systems, we recovered the planets discov¬ 
ered by the Kepler pipeline at high significance. For 
Kepler -42 and Kepler- 446, after removing the three sig¬ 
nals detected by the Kepler pipeline, there remained no 
significant peaks in the BLS, but for Kepler- 445, after re¬ 
moving the two known planet candidates, there remained 
a series of significant peaks spaced like harmonics of a 
transiting planet signal, shown in Figure [4j The highest 
peak in the BLS spectrum indicated a candidate period 
of 8.15275 days. Hereafter, we refer to this new planet 
candidate as Kepler- 445d, orbiting near a 5:3 resonance 
with Kepler- 445c, and we confirm all of the planet in Sec- 
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Fig. 4.— BLS periodogram of Kepler -445 with the two planets 
detected by the Kepler pipeline removed. There are a series of 
peaks with power greater than our threshold (horizontal blue line) 
corresponding to the true period (denoted by the vertical red hash 
mark) and harmonics of the new planet Kepler- 445d. 

tion [3J We also detected no significant photocenter shift 
in the centroid of the Kepler image of Kepler- 445 dur¬ 
ing the Kepler- 445d transit events. Taking the difference 
between the centroid of the image in and out of transit, 
we calculate a photocenter shift of 0.00054 ± 0.00033 
arcseconds. 

2.5.3. Transit Fits 

Given the large uncertainties on the physical radii of 
the planets orbiting Kepler-4,46, and the new detection of 
Kepler- 445d, we chose to fit the the transit light curves 
for Kepler-445 and Kepler-446 independently. We fitted 
transit light curves to the data using a modified version 
of the Transit Analysis Pac kage, or TAP, an IDL soft¬ 
ware package developed by IGazak et all (|2012f ) . TAP 
employs a Markov-Chain Monte Carlo algorithm within 
a Bayesian framework for determining transit parame¬ 
ters. In order to reduce the degeneracy between impact 
parameter and planet-to-star radius ratio, we modified 
TAP to impose a strict prior on t he stellar densi t ies f or 
Keple r- 445 and Kepler-446. As ISeager fe Mallen-Ornelasi 
(|2003ll showed, ingress/egress duration and full transit 
duration can be combined in such a way so as to de¬ 
termine the density of the host star, assuming knowl¬ 
edge of the planet’s orbital period, orbital eccentricity 
and longitude of periastron. Likewise, knowledge of the 
host star density, orbital eccentricity and longitude of 
periastron can be combined to constrain the relation¬ 
ship between ingress/egress duration and transit dura¬ 
tion. In the case of the planets orbiting Kepler- 445 and 
Kepler-446, ingress/egress duration is difficult to mea¬ 
sure due to the V-shaped nature of the light curves as 
discussed above, and imposing such a constraint is a 
powerful way to fit ac curate transit parameters. Indeed, 
iMuirhead et all (|2012f> used this technique to fit the tran¬ 
sit light curves for the Kepler- 42 system, although they 
fixed the stellar density rather than applying a prior. 

We assumed a stellar density prior based on the mea¬ 
sured masses and radii for Kepler- 445 and Kepler- 446 
from spectroscopy described in the previous section. Due 
to the use of evolutionary models, the mass and radius 
uncertainties are nearly 100% covariant for the stars, so 
we only use the uncertainty in mass when calculating the 


density prior. We assumed a Gaussian prior with a mean 
of 26.51 gm cm -3 and a standard deviation of 6.54 gm 
cm -3 for Kepler- 445, and a mean of 22.73 gm cm -3 and 
a standard deviation of 6.06 gm cm' 3 for Kepler- 446. 

We also assumed that all of the planets orbiting Kepler- 
445 and Kepler-446 have low eccentri city (e), an d fix 
the va lue to zero in our fit. Following IWu fe Goldreichl 
l|2002t , Eq. 1), and assuming reasonable planet densities 
with tidal Q' P values of 100 to 10000, all of the planets 
orbiting Kepler-445 and Kepler-446 have circularization 
timescales of less than 1 Gyr. Given the lack of quiescent 
Ha or u-band emission in either Kepler- 445 or Kepler- 
446, and Kepler- 446’s high galactic space motion, we are 
confident that the stars are older than 1 Gyr and likely 
older than 5 Gyr. We note, however, that planets can 
sustain eccentricities over long timescales via spin-orbit 
coupling, and that in multiple-planet systems non-zero 
eccentricity can be sustained via orbital resonances. The 
assumption of zero eccentricity also negates the impor¬ 
tance of the longitude of periastron parameter for the 
transit fit or for the imposed prior on stellar density. 

For limb-darkening, we fit a linear (ul) and quadratic 
(u2) coefficient to the transits. For each system, Kepler- 
445 and Kepler-446, we tied the limb-darkening param¬ 
eters across each of the transiting planet fits. We used 
Gaussian priors for the limb-darkenin g coe ffici ents, based 
on th e expected coefficients from IClaret fe Bloemenl 
(|2011f) using our measured stellar parameters and their 
uncertainties. For Kepler-445 we used a Gaussian prior 
with a mean linear limb-darkening parameter of 0.50 
with standard deviation of 0.17, and a quadratic limb- 
darkening parameter of 0.35 with a standard deviation 
of 0.13. For Kepler-446, we used a Gaussian prior with a 
mean linear limb-darkening parameter of 0.42 with stan¬ 
dard deviation of 0.12, and a quadratic limb-darkening 
parameter of 0.35 with a standard deviation of 0.11. 

Finally, we chose to fix the orbital periods and transit 
epochs (ffi) of the planets to the values listed in the NEA 
for all but Kepler-4456, for which we fit a separate tran¬ 
sit model to the Kepler light curve using a Levenberg- 
Marquardt minimization routine. We used the best fit 
period and transit epoch for the TAP fit. Allowing the 
periods of the planets to vary within TAP had no sig¬ 
nificant effect on the resulting transit parameters. Our 
final transit parameters for the six planets are listed in 
Tables [3] and [3 and we show the phase-folded Kepler 
light curves and fitted transit curves in Figure [5j We 
combine the stellar parameters with the transit parame¬ 
ters to determine the physical parameters for the planets, 
which are also listed in Table 0 and Figure [3] illustrates 
the planets with Kepler- 42 and the Galilean moons of 
Jupiter for comparison. We also calculate the incident 
flux on the planets as a fraction of the solar flux inci¬ 
dent on the Earth’s upper atmosphere (So = 1360 W 
m -2 ). The values indicate that all of the planets are 
likely too hot to be located within their host stars’ hab¬ 
it able_zones, using habitable-zone limits calculated by 
[Kopparapul (|2013lf . However, with an S/So of 2.24, one 
could argue that Kepler-4456 is near the habitable zone. 

We also include coarse estimates for the planet masses 
and expected semi-amplitude radial velocity signatures 
in Table [H using the recent em pirical ly-measured planet 
mass-radius relations of iMarcv et all (120141 1. All of the 
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TABLE 3 

Transit Parameters for Kepler- 445 


Parameter 

Kepler- 445b 

Kepler- 445c 

.R'ep/er-445d 

P (days) 3, 

t 0 - 2454833 (BJD) a 

Rp/R* 

b 

T (hours) 

Ml b 

M2 b 
e a 

2.984151 ± 0.000011 

133.1194 ± 0.0033 

0.0676 ± 0.0018 

0 001 + 0 - 203 

U.UU1 —0.001 

1.0304 ± 0.0180 
0.141 ± 0.080 
0.263 ± 0.071 

0 

4.871229 ± 0.000011 

133.6408 ± 0.0019 

0.1075 ± 0.0014 

0 000 + 0 ' 101 

U.UUU -0.000 

1.2287 ± 0.0154 
0.141 ± 0.080 
0.263 ± 0.071 

0 

8.15275 ± 0.00040 

3.7512226 ± 0.05 

0.0533 ± 0.0029 

0 011 +°- 470 

U.Uli -0.010 

1.3860 ± 0.0888 
0.141 ± 0.080 
0.263 ± 0.071 

0 

Inc (degrees) 0 
a/R* c 

89 74 +°- 18 
-0.28 

21.94 ± 0.30 

89.91 IHo 

30.21 ± 0.38 

89.61 ±° 0 2 2 l 

42.95 ± 0.58 

R P (R®y 

1.58 ± 0.23 

2.51 ± 0.36 

1.25 ± 0.19 

S/S 0 d 

8.57 ± 0.60 

4.52 ± 0.31 

2.24 ± 0.17 

M P (M®) e 

4 to 6 

8 to 9 

3 to 4 

K (m s 1 ) e 

6 to 8 

9 to 10 

2 to 4 


a Held fixed in fitting procedure. Period s and ephemerides fo r Kepler- 445c and Kepler- 445b 
are from the NASA Exoplanet Archive (jAkeson et al.lf2013T ). 

b Limb-darkening coeffic ients were tied between planets, and subject a prior based on 
I Claret fe Bloemenl J2011J) . 
c Calculated from parametrization. 

d So = 1360 W m —2 , the solar flux incident on Earth’s upper atmosphere. 
e Coarsely e s timat ed using the empirically-measured planet mass-radius relationships of 
IMarcv et al.l (|2014T ). 


TABLE 4 

Transit Parameters for Kepler -446 


Parameter 

Kepler- 446b 

Kepler- 446c 

Kepler- 446d 

P (days) 3 

to - 2454833 (BJD) 3 
Rp/R * 
b 

T (hours) 

Ml b 

M2 b 
e 3 

1.565409 ± 0.0000033 

132.9135 ± 0.0019 

0.0574 ± 0.0026 

0 601 + 0 - 096 

U.OUl -0.088 

0.6456 ± 0.0456 
0.447 ± 0.059 

0.353 ± 0.065 

0 

3.036179 ± 0.0000055 

134.069573 ± 0.000945 

0.0424 ± 0.0018 

0 025 +°- 529 
u.uzo -0.025 

0.9432 ± 0.0504 
0.446 ± 0.059 

0.353 ± 0.065 

0 

5.148921 ± 0.000022 
133.3196 ± 0.0042 
0.0519 ± 0.0022 

0 705 +°-°57 

u. IUO -0.066 

0.8832 ± 0.0480 
0.446 ± 0.059 

0.353 ± 0.065 

0 

Inc (degrees) 0 
a/R* c 

87.42 

14.20 ± 0.94 

88.97 

22.40 ± 1.36 

88.72 

31.60 ± 2.08 

Rp ( R<b) c 

1.50 ± 0.25 

1.11 ± 0.18 

1.35 ± 0.22 

S/S 0 d 

26.22 ± 4.70 

10.54 ± 2.03 

5.30 ± 0.83 

M P (lW®) e 

4 to 5 

2 to 4 

3 to 5 

K (m s" 1 ) 6 

6 to 8 

2 to 5 

3 to 5 


a H eld fixed in fittin g procedure. Periods and ephemerides are from the NASA Exoplanet Archive 
JAkeson et al.|[20130 . 

b Limb-darkening coeffi cients were tied between planets, and subject a prior based on 
I Claret Sz Bloemenl J2011J) . 

° Calculated from parametrization. 

d So = 1360 W m —2 , the solar flux incident on Earth’s upper atmosphere. 

e Coarsely estima ted using the empirically-measured planet mass-radius relationships of 
IMarcv et al.l J2014jj . 
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planets have anticipated radial velocity semi-amplitudes 
of over 1 m s -1 . However, being inid-M dwarf stars, 
the stars are relatively faint for current-generation 
visi ble-light pr ecisi on-radial-velocity spectrometers (e.g. 
IBottom et al.1 120131 1. In the future, the stars may 
be compelling targets for next-ge neratio n, infr ared 
preci si on-radial-velocity instruments dMa hadeva n et al.l 
2012t lH alverso n et all [2 014 iQuirrenbach et al.| 2012|; 


Barrick et al. I2012t lArtieau et al.1 20121: Thibault et al 


2012:lMicheau et ahll201 jlPares et al.ll2012l : lCrenD etaT 


2014 iGe et al.l 12014 1. Such measurements would 

prove useful for precisely measuring the planet masses 
and for placing constra ints on their atmospheres (e.g. 
iMiller- Ricci Kempton et ~ah I2012H and interior struc ¬ 
tures (e.g. iRoeers fe Seaeei 120101 : [Valencia et al.l 12013H . 
as w ell as measur in g any non-zero eccentricity (e.g. 
lAnglada-Escude et al.ll2013l f. 


3. FALSE POSITIVE ANALYSIS 

As with the vast majority of Kepler planet candidates, 
the transit signals detected in the light curves of Kepler- 
445 and Kepler- 446 are not amenable to dynamical con¬ 
firmation as bona fide planets, either by radial velocity or 
transit timing variation measurements. Confirming their 
planetary nature thus requires probabilistic validation; 
that is, demonstrating that the probability for them to 
be caused by a blended stellar eclipsing binary, or any 
other astrophysical false positive scenario, is very low. 

While broad arguments have demonstrated that only 
a small number of Kepler planet candidates are ex¬ 
pected to turn out to b e astrophysical fals e positives 
llMorton fe Johnsonl 120111 : iFressin et al.l l2013h , any spe¬ 
cific planet candidates of particular interest should be in¬ 
dividually investigated more thoroughly to ensure their 
planetary nature. To this e nd, we h ave employed the 
method described in detail bv lMorton ; 2012 i to calculate 
false positive probabilities (FPPs) for all six of the can¬ 
didates discussed in this paper. This method compares 
the shapes of the observed light curves to the shapes of 
simulated planet and astrophysical false positive scenar¬ 
ios in order to determine the relative likelihoods of the 
signals to be caused by the different scenarios. 

The basic assumption of this analysis is that the ob¬ 
served transit signal is roughly spatially coincident with 
the target star—that is, that the flux decrement in the 
aperture of the target star is not due, for example, to 
a different star several arcseconds separated from the 
target star. An important part of the Kepler vett ing 
pipeline is looking for such offsets dBrvson et al.ll2013h — 
only candidates that do not show significant offsets get 
promoted to “planet candidate” status. The results 
of these tests are summarized in the NASA Exoplanet 
Archive tables in the column titled “PRF A(9mq” and as¬ 
sociated uncertainty, which indicates the difference in po¬ 
sition between the Kepler pixel response function (PRF) 
fitted to the in- and out-of-transit data. For Kepler- 
446b, c and d and Kepler- 445b and d, all these offsets 
are significantly smaller than 1 arcsecond. For Kepler- 
445d, which is not in the Kepler catalogs, we perform 
our own in- and out-of-transit centroid analysis, finding 
no noticeable shift. Whil e we note that this analysis is 
not as detailed as the lBrvson et al.1 (120131) PRF analysis, 
it certainly indicates that there is no large offset of the 
candidate signal. 







Fig. 5. — Phase-folded light curves for the Kepler- 445 and Kepler- 
446 transiting planets with our best fit transit light curves. Raw 
data are shown as grey dots, binned data as blue circles, and the 
model as red lines. Note the V-shaped transits for Kepler- 446, 
which is due to the ~30 minute integration times of the Kepler 
data. We speculate that this contributed to the erroneous planet- 
to-star radius ratios returned by the Kepler pipeline. Also note 
that the scale of the y-axis changes for the Kepler- 445 and Kepler- 
446 plots. 
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Fig. 6.— Scaled illustration of the Kepler-42, Kepler-445 and Kepler -446 exoplanetary systems, with Jupiter and the Galilean moons 
shown for reference. The colors of the stars match what is expected given the visible-light spectra shown in Figure[2] and the limb-darkening 
applied to the images of the stars matches the measured limb-darkening from the transit curves. The dramatic difference in the color and 
limb-darkening between Kepler -445 and Kepler-446/Kepler-42 is likely due to the markedly different metallicities, though we stress these 
are all M4 dwarf stars by spectral classification. Interestingly, the higher metallicity of Kepler -445 results in a redder r-J color, but a bluer 
appearance to the human eye compared to metal-poor stars, due to the peculiar wavelengths of deep molecular features within the visible 
spectrum (400 to 700 nm). 


While the pixel-level analysis of these candidates lo¬ 
calizes them to be at least within 1 arcsecond of the 
target stars’ locations, limiting the possible on-sky area 
in which could exist a potential false-positive blend, the 
constraints from the AO and aperture masking sensitiv¬ 
ity curves presented in Section 12.11 further restrict the 
parameter space of possible fals e positiv es. Incorporat¬ 
ing all these constraints into the lMortonl (l2012f) analysis, 
we obtain FPPs for all candidates of < 10 -4 , except for 
Kepler-AAbd and Kepler- 446d, for which we obtain FPPs 
of approximately 1 in 500 and 1 in 200, respectively. 
For both Kepler- 445d and Kepler-AAdd the most prob¬ 
able false positive scenario is an eclipsing binary within 
a hierarchical triple stellar system. 

We note, however, that these FPP calculations do not 
account for the fact that the candidates are observed to 
be in multiple-transiting systems. As it is now known 
that many multiple planet systems have low mutual in - 
clinations (iFabrvckv et all l2012t iFang fc Margotl 120121 ) , 
the presence of one transiting planet means that addi¬ 
tional planets have a significantly higher probability to 
transit than if their orbital inclinations were randomly 
distributed. Accounting for this effect gives a “multi¬ 
plicity boost” to the probability of the planet scenario 
of about a factor of ~5-10 (depending on the exact as¬ 
sumptions) , thus lowering the FPPs of both Kepler- 445d 
and Kepler-AAdd to below 0.1%— all six of the planet 
candidates discussed in this paper are probabilistically 
validated. 

However, there is another scenario, while not strictly 
an astrophysical false positive, that we cannot completely 
rule out. That is, the AO and aperture masking ob¬ 
servations presented here do not completely exclude the 
presence of stellar companions around Kepler- 445 and 
Kepler-446. This can be understood by noting that the 
inner working angle of 20 mas achieved by the aper¬ 
ture masking observations corresponds to a projected dis¬ 
tance of 2 AU at lOOpc, and that a snapshot of a binary 
star system will generally not be at maximum separa¬ 
tion, due to orbital position and inclination. Using a 
Monte Carlo simulation of potential stellar companions 
distributed according to the stellar binary period and 


eccentricity distribution of lRaghavan et all (120101) (with 
randomized orbital elements), we find that about 25% 
of potential stellar companions would be undetected by 
these high-resolution imaging observations. However, a 
close stellar companion would be detrimental to the ex¬ 
istence of the observed planets, so we place another con¬ 
straint requiring putative stellar companions to have pe¬ 
riods longer than 200 days; this removes another ^15% 
of potential stellar companions. Thus, starting with an 
assumption that ab out 40% of M-dwarf stars are in stel¬ 
lar b i nary systems (iFischer fo Marcvl 1 1 9!)2t iClark et al.1 
120121 : iDuchene fe Krausll2013lh then our observations al¬ 
low about a 4% chance for either of Kepler-445 or Kepler- 
446 to be binary s ystems. _ 

As discussed in [Muirhead et al.i (120121 b there would 
be two potential consequences of the presence of such 
a stellar companion in either of these systems. First, 
the planet radii would be larger than we estimate due 
to the unaccounted-for dilution (by a maximum of \/2 
if the planets are around the primary star, and by a 
larger factor if they are around a fainter secondary star). 
Second, there is a possibility that the transiting plan¬ 
ets could be distributed around both stars in the system, 
rather than all being around the same, single star. Such a 
“split-multi ” co nfiguration, while predicted to be rare by 
iFabrvckvl (120111b does undou btedly exist in the arc hetype 
KOI-284/Kepler-132 system dLissauer et al.ll2014lb a vi¬ 
sual binary with two of its three candidates having orbital 
periods of 6.18 and 6.42 days. However, given the low 
probability for either of these two systems to be binaries, 
and the even lower probability that the planets would be 
split between the two components, we neglect this sce¬ 
nario for the remainder of this paper, although we do 
acknowledge its possibility at the ^^-probability level. 

4. OCCURRENCE OF COMPACT MULTIPLES ORBITING 
MID-M DWARF STARS 

With three confirmed compact multiple systems or¬ 
biting rnid-M dwarf stars-Jbepfer-42, Kepler-445 and 
Kepler- 446-we can estimate the occurrence of compact 
multiples like these around mid-M dwarf stars in gen¬ 
eral: defined here to be a nearly coplanar system of 2 
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Fig. 7.— Top : Density contours for KIC J — H vs. r — J of 
all objects observed by Kepler for at least one quarter from Q1 
to Q 16, with spectrosc opically confirmed dwarf and giant stars 
from [Mann et all (|2012l . , green and red circles, respectively). Bot¬ 
tom: Same, but with J-K , which provides better dwarf-giant dis¬ 
crimination with no confirmed giant stars contaminating the dwarf 
sequence. To calculate the occurrence of compact multiples or¬ 
biting mid-M dwarf stars, we isolate mid-M dwarf stars using a 
the color cuts indicated in the lower plot. Kepler- 42, Kepler- 445 
and Kepler- 446 are shown ( blue triangles ), using our independent 
r magnitude measurements. We include GJ 1214, as it is the only 
other mi d-M dwarf know n to host a confirmed transiting planet 
ICharbonneau et al.l (120091) . The lower positions of Kepler- 446 and 
Kepler-42 compared to Kepler- 445 and GJ 1214 in J-H are consis¬ 
tent with their relative metallicities as measured by infr ared spec¬ 
troscopy dRoias- Ayala et al.II2013 ; I Muir head et aLll2014j) . We note 
that many more mid-M dwarf stars were added in quarter 17, but 
this quarter was cut short and was not included in the current exo¬ 
planet search by the Kepler pipeline. Therefore, we did not include 
them in the occurrence calculation. 


or more planets all orbiting with periods of less than 10 
days. First, we downloaded and combined the lists of all 
targets observed by Kepler for at least one quarter from 
quarter 1 to 16 from MAST0. The files include measure¬ 
ments o f the combined differenti al photometric precision 
1CDPP. [Christiansen et ahll2012f) for each target for each 
quarter, calculated over 3, 6 and 12 hour durations. 

To calculate the occurrence of compact multiples orbit¬ 
ing mid-M dwarf stars, we u se two approac h es. The first 
follows the methodology of IHoward et al.l (|2012t ). who 


http://archive.stsci.edu/pub/kepler/catalogs/ 


calculated planet occurrence as a function of planet ra¬ 
dius and orbital period using the Kepler planet candi¬ 
date discoveries. We considered the number of mid-M 
dwarf stars around which Kepler and the Kepler pipeline 
could have found compact multiples like the ones orbiting 
Kepler-42, Kepler-445 or Kepler- 446. We exclude Kepler- 
445d from our calculation, as that was not discovered by 
the Kepler pipeline. We isolated mid-M-dwarf stars from 
the Kepler targets applying two color cuts using magni¬ 
tudes from the KIC: r— J > 3.2 to isolate red objects, and 
J — K < 0.0555 (r-J) + 0.7622 to remove evolved stars 
(giant and highly-reddened sub-giant stars). At first we 
considered using J — H to remove giant sta rs; h owever, 
we found that several giant stars identified in lMann et al.1 
(l2012f) were not excluded using J—H, but were excluded 
using J—K, so we chose the latter. We also excluded tar¬ 
gets that were only observed in the very first quarter of 
Kepler (quarter 0), which was only 30 days and has lower 
intrinsic photometric precision than subsequent quarters, 
and targets that were added in quarter 17, since that has 
not yet been searched by the Kepler pipeline and was cut 
short due to a spacecraft malfunction. The color cuts re¬ 
sulted in 509 mid-M dwarf stars with at least one full 
quarter of Kepler observations. Figure [7] illustrates the 
color-color cuts, and the position of Kepler-42, Kepler- 
445 and Kepler- 446. 

Of the 509 mid-M dwarf stars, Kepler-42, Kepler-445 
and Kepler -446 are the only stars with multi-planet- 
candidate systems listed on the NASA Exoplanet Archive 
at the time of this study. In addition to these three 
stars, KOI-2862, KOI-4290, KOI-3855, KOI-3138 and 
KOI-5327 also host single transiting planet-candidates 
and fall within the defined color cuts. We do not con¬ 
sider the single-transiting-planet hosts in this occurrence 
calculation. 

We then calculated the number of mid-M dwarf stars 
with the signal-to-noise ratio necessary to discover each 
of the planets orbiting Kepler- 42, Kepler -445 or Kepler- 
446. To calculate the signal-to-noise ratio (S/N) for an 
arbitrary transiting planet orbitin g a partic ular Keple r 
target, we used Equation A1 from lFressin et all (I2013T ). 
repeated here with a slight modification: 

s / n =c5^# 

where <5 is the transit depth, CDPP e g is the combined 
differential photometric precision over the duration of 
the transit, fobs is the total Kepler observing time used 
to search for the planet by the Kepler pipeline, and P 
is the planet orbital period. To determine CDPP e ff, 
we first averaged the 3, 6 and 12-hour CDPP values 
for each mid-M dwarf target over all quarters, then 
fit a function to the CDPP versus time with the form 
CDPP e ff = CDPPo/\/tdur > where fdur is the duration 
of a transit. We then calculated the signal-to-noise that 
would be achieved on each of the planets orbiting Kepler- 
42, Kepler-445 (excluding Kepler- 445d) and Kepler- 446, 
if they were orbiting every single mid-M dwarf observed 
by Kepler with the same transit durations. 

We adopt a value of 10 as the minimum signal-to- 
noise r equired to detect a planet, following IHoward et al.1 
(|2012T) . This is a reasonable value, considering the signal- 
to-noise needed to detect the Kepler- 42, Kepler-445 
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TABLE 5 

Planet S/N and Transit Probabilities 


Planet 

Transit S/N 
for disposition 

Kepler Mid-Ms 
where S/N 
would be >10 

Transit 

probability 

Kepler -42 b 

67.9 

353 

0.0826 

Kepler-42 c 

94.3 

392 

0.1587 

Kepler-42 d 

35.2 

261 

0.0621 

Kepler- 445c 

26.3 

468 

0.0331 

Kepler- 445b 

12.8 

431 

0.0456 

Kepler- 446b 

24.5 

412 

0.0704 

Kepler- 446d 

10.2 

300 

0.0316 

Kepler- 446c 

14.5 

301 

0.0446 


and Kepler -446 transit signals and give them a planet- 
candidates disposition. For the three planets orbiting 
Kepler-A2 , they were detected by the Kep ler p ipeline us¬ 
ing only the first quarter of data (jBorucki et al.lf20TT ). 
For Kepler-AA5 and Kepler- 446, we combine the “disposi¬ 
tion provenance” listed on the NASA Exoplanet Archive 
with the number of observed quarters within the prove¬ 
nance boundaries to determine the number of quarters 
required to detect their orbiting planets. Following this, 
Kepler- 445c and Kepler-AAbb were detected by the Ke¬ 
pler pipeline using two quarters of data (6 and 8). The 
three planets orbiting Kepler -446 were detected by the 
Kepler pipeline using 1 quarter of data (7) for 0.01 and 
0.02, and two quarters (7 and 12) for 0.03. Table [5] 
indicates the resulting S/N values for the eight plan¬ 
ets achieved when they were detected, all of which are 
greater than 10. We also list the number of Kepler mid- 
M dwarf stars for which the S/N of the respective planet 
would be higher than 10, if it were transiting all of the 
Kepler mid-M dwarf stars. Lastly, we list the individual 
transit probabilities (i?*/a). 

When calculating the number of mid-M dwarf targets 
for which Kepler- 42, Kepler -445 or Kepler- 446-like sys¬ 
tems could be detected, we make two assumptions. First, 
we assume the mid-M dwarf targets have similar radii to 
Kepler- 42, Kepler -445 and Kepler- 446, all roughly 0.20 
Rq. The true stellar radii of the mid-M dwarf targets 
affects the signal-to-noise of the Kepler-A2 , Kepler-AAb 
or Kepler-AAQ-Y\ke systems if there were transiting other 
stars. The assumption of similar radii across all the mid- 
M dwarf targets is a reasonable one considering that 
Kepler- 42, Kepler-AAb and Kepler -446 themselves have 
similar radii and nearly subtend the full color-color space 
defined in the mid-M dwarf target cut (see Figure 0. 
Even so, without spectra for the full sample of mid-M 
dwarf stars we must make this assumption. Second, we 
assume that the impact parameter has a negligible im¬ 
pact on the signal-to-noise of the planet transits if they 
were transiting other mid-M dwarf targets. In fact, an es¬ 
pecially large impact parameter will reduce the signal-to- 
noise of a transit event since the duration will be shorter 
and the transit could even be grazing. 

To calculate the oc currence o f com pact multiples, we 
use Equation 2 from iHoward et al.1 (|2012f ) , who calcu¬ 
lated the frequency of planets for a specific range of or¬ 
bital periods and planet-sizes, but modified for compact 


multiple systems: 


fT'sys 

/cm = ^2 

3 


1 /pJ 

sys,j 


( 2 ) 


Whereas IHoward et al.l (|2Q12h calculated the number 
of planets per star for specific ranges of orbital period 
and planet radius, we, on the other hand, want to cal¬ 
culate the fraction of mid-M dwarfs that host compact 
multiple systems, which we call /cm- Instead of index¬ 
ing each planet detection in the summation, instead we 
index each system: j refers to Kepler- 42, Kepler-AAb and 
Kepler-AAQ, successively, pj refers to the transit probabil¬ 
ity of the whole compact multiple (corresponding to the 
transit probability of the outermost system), and n sys j 
is the number of stars around which the system could be 
detected. Equation [2] implicitly defines a compact multi¬ 
ple as a system as being like Kepler-A2, Kepler-AAb and 
Kepler-AAQ: 2 or more planets all orbiting with periods 
of less than 10 days. 

Also implicit in this formalism is the assumption that 
compact multiple systems are co-planar: if the outer¬ 
most planet is detected, the inner planets would also be 
detected, presuming their signal-to-noise is high enough. 
If, in fact, compact multiple systems are not coplanar, 
then /cm represents a lower limit to their frequency 
around mid-M dwarf stars, as the single transiting- 
planet-candidate hosts metioned above could harbor 
many more interior planets that do not transit. 

We calculated the uncertainty in /c m fo l lowin g the 
same procedure as described in lHoward et al.l (120121) : We 
calculated the binomial probability distribution of draw¬ 
ing compact multiple planets from n v \/f stars, and use 
the closest values to ltx above and below /cm of the bino¬ 
mial distribution as the uncertainties (~ 15% and ~ 85% 
of the cumulative binomial distribution). 

Performing this calculation we arrive at a compact- 
multiple occurrence of 23 > 3 %. The uncertainties cor¬ 
respond to drawing 5 and 2 compact multiple systems 
from 14 stars, based on the binomial distribution. We 
note that we include mid-M dwarf stars of all metallici- 
ties observed by Kepler in this calculation, which roughly 
matches t he metallicit y distribution of the solar neigh¬ 
borhood ( Mann et al J i2013b ). and that Kepler-AAb has 
significantly super-solar metallicity ([Fe/H]=+0.25) and 
Kepler-AAQ and Kepler-A2 both have sub-solar metallic¬ 
ity ([Fe/H]<-0.30). We conclude that roughly one-fifth 
to one-quarter of mid-M dwarf stars host compact mul¬ 
tiple systems, for the full spectrum of metallicities in the 
solar neighborhood. 


4.1. A Monte Carlo Approach 

As a consistency check, we also performed the occur¬ 
rence calculation using a Monte Carlo approach, taking 
into account the effect of transit duration on SNR, as 
well as the probability of a Aepler-pipeline detection for a 
given SNR. We assigned each of the Kepler mid-M dwarf 
targets with at least one-quarter of data a random in¬ 
clination, then calculated the SNR of the transit signals 
if the Kepler- 42, Kepler-AA5 or Kepler-AA6 planets were 
transiting at the random inclina tion. We then use d the 
detectability ramp described in iFressin et al.l (120131) to 
randomly determine whether each system would be de- 
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Fig. 8. — Distribution of /cm calculated for Sun-like stars, early 
M dwarf stars and mid-M dwarf stars, subject to arbitrary normal¬ 
izations for clarity. For mid-M dwarf stars, we simulated Kepler- 
42-, Kepler- 445- and Kepler-446-\ike planets orbiting mid-M dwarf 
targets observed by Kepler for at least one full quarter, with ran¬ 
dom inclinations and subject to a randomized detectability, /cm 
was calculated using Equation [3] The distribution has a median of 
26%, a mode 21% (with 1% bins), and asymmetric uncertainties 
of %, corresponding to the equivalent of +/- la, which is con¬ 
sistent with the value determined by Equation \2\ We performed 
identical calculations for Sun-like stars and early M dwarf stars, 
both of which show lower values of /cm compared to mid-M dwarf 
stars, but the differences are less than 3a. 

tected by the Kepler pipeline. This process results in 
the number of Kepler-42-, Kepler-445- and Kepler- 446- 
like systems that should have been detected if they were 
orbiting all mid-M dwarf stars. Since one Kepler- 42, one 
Kepler- 445 and one Kepler- 446 system was detected, we 
use the following summation to determine /cm: 

/cm =- 1 - 1 - (3) 

^Kepler— 42 ^Kepler— 445 ^Kepler— 446 

where n refers to the number of systems that would have 
been detected if they orbited all mid-M dwarf stars. We 
repeated the simulation 10000 times, and the resulting 
distribution of /cm is shown in Figure [8] The distribu¬ 
tion has a median of 22%, a mode 21%, and asymmetric 
uncertainties of /g %, corresponding to the equivalent of 
+/- lcr, which is consistent with the value determined by 
the first method. 

We note that planets too small to be detected are not 
included in this occurrence calculation. When including 
planets smaller than the Kepler- 42 planets within the 
definition of “compact multiples,” this result becomes a 
lower limit on the occurrence of compact multiples orbit¬ 
ing mid-M dwarfs. 

4.2. Comparison to Sun-like and Early M Dwarf Stars 

We compared this result to the occurrence of compact 
multiples orbiting Sun-like stars and early M dwarfs. To 
create a sa mple o f Sun-like stars , we applied the crite¬ 
ria used by iPetigura et al.l (|2013l) to the latest catalog 
of stellar paramete rs for the Kepler targets, compiled by 
iHuber et al.l (l2014f) . We selected all Kepler targets with 


effective temperatures between 4100 and 6100 K, with 
surface gravities log(g) between 4.0 and 4.9, and with 
Kepler magnitudes between 10.0 and 15.0. We cross- 
referenced this list of Sun-like Kepler targets with the 
NASA Exoplanet Archive list of Kepler planet candi¬ 
dates and identified 93 systems meeting the criteria for 
a compact multiple: at least two transiting planet can¬ 
didates both orbiting with periods of less 10 days. The 
systems contain transiting planet candidates with phys¬ 
ical radii ranging from 0.45 1?® to 7.0 1?©. With a list 
of compact multiple systems found in a sample of Sun¬ 
like target stars, we repeated the calculation described 
above. However, we included a term in the SNR cal¬ 
culation that corrects for the differences in stellar radii 
between the compact multiple ho sts and target star s, us¬ 
ing the stellar radii reported in iHuber et all (i2014f ). We 
calculated a compact multiple occurrence of 8.8/ 0 ;g% for 
Sun-like stars, 2.7 ct lower than the occurrence for mid-M 
dwarf stars. 

As illustrated in Figure [3 early M dwarfs and evolved 
stars overlap in the r — J vs. J — K color-color diagram. 
Instead of using colors to isolate early M dwarf Kepler 
targets, instea d we used t h e list of M dwarf targets from 
(Dressing fc Charbonneaui 1 20131 . their Table 4), exclud¬ 
ing any M dwarfs with r — J > 3.2. We cross-referenced 
the list to the NASA Exoplanet Archive and found 13 
systems with transiting compact multiple candidates: 
KOI-248, KOI-251, KOI-571, KOI-898, KOI-899, KOI- 
936, KOI-952, KOI-1078, KOI-1681, KOI-1843, KOI- 
1867, KOI-2036 and KOI-2793. As with the Sun-like 
star calculation, we included a stellar radius correction 
due to the wide range of radii that early M dwarf stars 
can have (~0.3 to 0.6 Rq), taking the stellar radii from 
iDressing fc Charbonneaui (12013! ). The compact multiple 
planet candidates orbiting early-M dwarfs range in radii 
from 0.73 to 2.69 1?©, similar to the sizes of compact 
multiple planets orbiting mid-M dwarfs. Repeating the 
calculations described above but for early M dwarfs, we 
find a compact multiple occurrence rate of 15.9/ 1 ; 5 %, 
not significantly different from the occurrence from mid- 
M dwarf stars. 

Figure [8] depicts the resulting Monte Carlo distribu¬ 
tions of /cm for Sun-like stars, early-M dwarf stars and 
mid-M dwarf stars. Although the measured increase in 
compact multiple occurrence with later spectral type is 
not statistically significant, it may simply be a conse¬ 
quence of the increase in overall pla net occu r rence aro und 
M dwarf stars first calculated bv [Howard et al.l ( 2012 :). 

5. DISCUSSION 

Interestingly, Kepler- 445c is somewhat similar in na¬ 
ture to GJ 1214b: both are super-Eartlis or mini- 
Neptunes orbiting metal-rich mid-M dwarf stars. Kepler- 
445c o r bits w ith a longer period: 4.87 versus 1.58 days 
(jCarter et al.ll20lil) , and is subject to much lower stel¬ 
lar irradiation. It would be illuminating to study the 
atmosphere of this planet via transit transm issio n spec- 


troscopv, as has been done for GJ 1214b (Bean et al. 

2010 

: Croll et al J 

2011 Desert et al. 20111: Berta et al. 

2012 

JFraine et al.l 

2013; Kreidberg et al.ll2014D. if not for 


the star’s faintness: with a Ks band magnitude of 12.61, 
transit-transmission spectroscopy will be challenging. 


Both the Kepler-446 and Kepler-42 systems present 
cautionary tales for future transit surveys such as 
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NASA’ s Tra nsit Exoplanet Survey Satellite (TESS, 
iRickeif I2014T) . The Kepler pipeline assigned unphysi- 
cally large planet-to-star radius ratios with significant 
uncertainties for the planets in both systems. In the fu¬ 
ture, systems like these may have accurate astrometric 
parallax measurements from ESA’s Gaia Mission. Com¬ 
bining the Gaia absolute magnitudes with mass-radius- 
luminosity relations would enable accurate stellar density 
priors to be employed during the pipeline transit fitting 
algorithm. We suggest such a procedure, as TESS is ex¬ 
pected to observe significantly more M dwarf stars than 
Kepler. 

Following the same calculation in the introduction per¬ 
formed on Kepler-42, but for Kepler-445 and Kepler- 446, 
we calculate that their “1% of M*” protoplanetary disks 
consisted of 16.8 and 8.0 M© of metals, respectively. In 
the case of Kepler- 446, where the three planets have radii 
less than or equal to 1.5 i?©, we speculate that the three 
planets have primarily rocky compositio ns based on the 
rocky/non-rocky threshold calculated bv iRoeersl (12014) ; 
however, we note that the planet radii are very near the 
rocky/non-rocky threshold. Assuming all are rocky, we 
calculate that the three planets constitute 6.5 M© of met¬ 
als, meaning that ~80% of the Kepler- 446’s disk metals 
went into these three planets. We repeat that this calcu¬ 
lation relies on simplistic assumptions, but could indicate 
that rocky planet formation is efficient around low-mass 
stars with compact multiple exoplanets. 

The same calculation for Kepler- 445 is trickier. At 2.51 
i?©, Kepler- 445c is very likely to have a thick gaseous en¬ 
velope. However, at 1.58 and 1.25 i?©, Kepler- 445b and 
Kepler- 445d could be rocky, with masses of 3 .9 and 1.6 
M©, when interpolating onto the iFortnev et al.l (120071) 
relations for primarily rocky composition. That leads 
to 30% of the metal content of Kepler- 445’s protoplan- 
ety disk going into its orbiting planets, not accounting 
for the core of Kepler- 445c, which is presumably more 
massive than Kepler- 445b in order to have accreted gas 
during the protoplanetary disk phase. Assuming Kepler- 
445c has a rocky core as massive as Kepler- 445b, this 
leads to 56% of the metal content in these three planets 
alone. 

These results can be interpreted as one or a combina¬ 
tion of the following: (1) rocky planet formation is highly 
efficient around mid-M dwarf stars with compact multi¬ 
ple planetary systems, (2) disk mass fractions around M 
dwarf stars are significantly higher than 1% of the host 
star mass, and/or (3) the planets presented in this work 
are not in fact rocky, with significant mass in the form of 
hydrogen and helium. As stated earlier, next-generation 
near-infrared precise-radial-velocity spectrometers may 
be able to confirm the rocky nature of these planets. If, 
in fact, rocky planet formation is efficient for compact 
multiple systems, we would not expect to find signifi¬ 
cantly more planets beyond the orbits of Kepler- 42 d, 
Kepler- 445c or Kepler- 446c, simply because there is not 
enough rocky material in the protoplanetary disk to cre¬ 
ate significantly more planets. Therefore, we anticipate 
a dearth of outer, long-period planets with radii > 0.8 
i?© orbiting mid-M dwarf stars with compact multiple 
systems. 

The planet-metallicity correlation, wherein higher- 
metallicity stars are more likely to host gaseous 
planets, has been shown to be true for sun-like 


stars as well as early M dw ar f stars dSantos et alJ 
20011 iFischer fc Valentil 120051 I Johnson et al.l 120101 : 
Roias-Avala et al .1120 1 011 . However, it has not been thor¬ 
oughly investigated for mid- M dwarf stars due to the 
challenges in finding planets around them. It is there¬ 
fore intriguing that Kepler-445 ([Fe/H]=+0.27±0.13 
Muirhead et all l20H and GJ 1214 ([Fe/H]=+0.20 
Roias-Avala et al.ll2012f) are the only mid-M dwarf stars 
with transiting planets larger than 2i?©, and both 
stars are metal-rich compared to the sol ar neighbor¬ 
hood l[Fc/H]=-0.05 iGaidos fc Mannl 120 lH ). Based on 
this alone, we speculate that the planet-metallicity cor¬ 
relation may in fact extend down to mid-M dwarf stars 
for large planets (> 27?©). 
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